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Introduction: Understanding how the distribution 
of boulders on the lunar surface changes over time is 
key to understanding small-scale erosion processes and 
the rate at which rocks become regolith. Boulders de- 
grade over time, primarily as a result of micrometeorite 
bombardment [1,2], so their residence time at the sur- 
face can inform the rate at which rocks become rego- 
lith or become buried within regolith. Because of the 
gradual degradation of exposed boulders, we expect 
that the boulder population around an impact crater 
will decrease as crater age increases. Boulder distribu- 
tions around craters of varying ages are needed to un- 
derstand regolith production rates, and Lunar Recon- 
naissance Orbiter Camera (LROC) Narrow Angle 
Camera (NAC) images [3] provide one of the best 
tools for conducting these studies. Using NAC images 
to assess how the distribution of boulders varies as a 
function of crater age provides key constraints for 
boulder erosion processes. 

Boulders also represent a potential hazard that must 
be addressed in the planning of future lunar landings. 
A boulder under a landing leg can contribute to deck 
tilt, and boulders can damage spacecraft during land- 
ing. Using orbital data to characterize boulder popula- 
tions at locations where landers have safely touched 
down (Apollo, Luna, Surveyor, Chang’e-3) provides 
validation for landed mission hazard avoidance plan- 
ning. Additionally, counting boulders at legacy landing 
sites is useful because: 1) LROC has extensive cover- 
age of these sites at high resolutions (~0.5 m/pixel). 2) 
Returned samples from craters at these sites have been 
radiometrically dated, allowing assessment of how 
boulder distributions vary as a function of crater age. 
3) Surface photos at these sites can be used to correlate 
with remote sensing measurements. 

Most boulder populations occur in association with 
steep slopes as a result of mass wasting or around 
young, fresh craters that are large enough to have ex- 
cavated bedrock [4-7]. Here we use NAC images to 
analyze boulder distributions around three craters of 
varying sizes and ages: Cone crater (CC; 26 Ma [8], 
340 m diameter) at the Apollo 14 site (Fig. 1), North 
Ray crater (NR; 50 Ma [8], 950 m) at Apollo 16, and 
Surveyor Crater (SC; 200 Ma [8], 200 m) at Apollo 12. 

Methods: We use CraterTools [9] in ArcMap to 
visually identify and estimate the size of boulders. In 
NAC images, boulders are positive relief features and 
appear as bright, sun-facing pixels adjacent to dark, 
generally elongated shadows in low-sun images. Boul- 


der sizes are recorded in terms of a circular diameter to 
capture the longest dimension. Using NAC images 
with a 0.5 m/pixel resolution, the smallest boulders 
that can be identified with confidence are ~1 m (>3 
pixels including the shadow). We determine the dis- 
tance of each boulder from the center of the crater us- 
ing the haversine formula. 

Boulder Distributions: The distributions we inves- 
tigate are boulder size-range distributions, size- 
frequency distributions, and frequency-range distribu- 
tions. These distributions are important for understand- 
ing how far craters of varying sizes and ages distribute 
boulders. We determine boulder distributions at in- 
creasing distances (in units of crater radii) to find how 
the frequency of boulders varies as a function of dis- 
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Fig. 1: Top: Boulder counts centered at Cone crater. Colored 
circles and numbers indicate distances from the rim (black 
circle) in crater radii. Bottom: Size-range distribution of 
boulders around Cone. Larger boulders lie closer to the rim 
and few boulders are seen beyond 8 crater radii. 


tance from the crater rim. We omit boulders inside the 
rim, because steep slopes inside rims tend to refresh 
the rock population as crater degradation proceeds. 

Results: Our boulder distributions show that the 
largest boulders within a count area occur closer to the 
crater rim (within ~2-4 crater radii) and smaller boul- 
ders occur at all distances (Fig. 1). The largest boul- 
ders seen are ~8 m at Cone, ~22 m at North Ray, and 
~4 m at Surveyor. The quantity (areal density) of boul- 
ders decreases with increasing distance from the crater 
rim. The presence of boulders drops off around 8 crater 
radii for both CC and NR. Fewer boulders are seen 
around SC, and the size- and range-frequency distribu- 
tions are irregular compared to NR and CC. 

Implications: Because boulders degrade over time, 
in principle, more boulders should be present around 
younger craters. Previous studies have shown a strong 
relationship between crater age and ejecta rockiness 
[10,11], with m-scale boulder survival rates at craters 
<1 km in diameter on the order of ~150-300 Ma [10]. 
Our counts agree with this conclusion so far; fewer 
boulders are present at the older Surveyor crater (200 
Ma) than at the younger CC and NR craters. CC and 
NR are young relative to the estimated time required to 
break down boulders [10,11], so most of their boulders 
are still present on the surface. Cumulative rock abun- 
dances around craters of varying ages will help con- 
strain the boulder survival rates more accurately. 

Size is also important — larger craters tend to exca- 
vate more boulders than smaller craters, and the size of 
the largest boulder ejected is likely related to the size 
of the crater [12-15]. We see both a higher number and 
larger boulders present around the larger NR crater 
than around CC and SC. Power-law fits to the size- 
range and -frequency distributions, of the form y=ax”, 
for craters of varying size will provide enough a and b 
values to allow us to derive equations to predict these 
values as a function of crater size. These derived equa- 
tions can then be used to estimate the maximum boul- 
der size expected at a crater of a given size. 

The presence of fewer boulders and the irregular 
frequency distributions at SC are likely owing to its 
location in the vicinity of several other fresh craters 
with their own boulder distributions, as well as to its 
smaller size and older age. The boulders that are pre- 
sent around SC are much smaller than those seen at the 
two younger craters, owing to its smaller diameter. 
Counts at more craters with varying sizes and ages 
(and far from other fresh craters) will allow us to better 
constrain boulder survival rates and test quantitative 
models set forth by previous studies [i.e., 11]. 

Surface Photography: The location and size of 
boulders imaged and sampled by Apollo astronauts are 
very well documented, enabling comparison with and 
verification of our NAC measurements. We have used 
Apollo 14 and Apollo 16 surface photographs and 


boulder measurements [16,17] to locate and validate 
size measurements of boulders. 

Diviner Rock Abundance: NAC boulder distribu- 
tions can also be used to validate Diviner rock abun- 
dance (DRA) data, which measures the areal density of 
the surface covered in boulders [4]. We calculate the 
cumulative areal fraction (CAF) of the surface covered 
in boulders at Cone and North Ray to compare with the 
RA values presented in [4], and find that our LROC 
counts closely predict the same CAF as DRA values. 

Landing Site Safety: Finally, the results from boul- 
der counts on ejecta blankets will provide constraints 
on how far craters of varying sizes distribute boulders. 
Predicting boulder size distributions as a function of 
distance from a crater is particularly useful in inform- 
ing potential boulder hazards for future missions. 
Boulder distributions can be used to estimate and au- 
tomate counts for future missions using the power-law 
relationships derived from size-range distributions. 
Coupling LROC data with DRA also allows us to ex- 
trapolate boulder count trends to sub-meter boulder 
populations that are unresolvable with NAC images 
and that may also pose a landing hazard [18]. 

Future Work: Ongoing work includes conducting 
more counts of boulder populations around craters of 
varying ages and testing key parameters such as crater 
size, degradation state, terrain type (mare vs. non- 
mare), and regolith thickness on boulder populations. 
This work will allow us to more comprehensively as- 
sess how boulder distributions change as a function of 
time, thereby informing the median survival time of 
boulders and the Moon’s regolith production rates. 
When coupled with counts at other spacecraft landing 
sites and verified using surface photography, these 
results can also inform boulder populations at varying 
distances from craters and aid in establishing safe land- 
ing zones for future missions [13,14]. 
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